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ABSTRACT 

We introduce the tracer mass estimator. This is a new and simple way to esti- 
mate the enclosed mass from the projected positions and line of sight velocities of 
a tracer population (such as globular clusters, halo stars and planetary nebulae). 
Like the projected mass estimator, it works by averaging (projected distance) x 
(radial velocity)^/^ over the sample. However, it applies to the commonest case 
of all, when the tracer population does not follow the overall dark matter density. 
The method is verified against simulated datasets drawn from Monte Carlo re- 
alisations of exact solutions of the collisionless Boltzmann equation and applied 
to the recent M31 globular cluster data set of Perrett et al. (2002), as well as to 
M31's satellite galaxies. 

Subject headings: galaxies: general - galaxies: haloes - galaxies: kinematics and 
dynamics -galaxies: individual: M31 - dark matter 



1. Introduction 

Simple estimators based on the virial theorem have long been used to ascertain the 
masses of galaxy groups from radial velocities and projected positions of the members (e.g. 
Zwicky 1933). Nowadays, they have been largely superceded by the projected mass estimator, 
which is based on averaging (projected distance) x (radial velocitjY/G over the sample 
(Bahcall & Tremaine 1981; Heisler, Tremaine & Bahcall 1985). In the form in which it is 
usually encountered, a key assumption underlying the estimator is that the member galaxies 
track the mass of the group. The problem of estimating the masses of elliptical galaxies 



- 2 - 



and the haloes of spiral galaxies is different. Here, the datasets available are the radial 
velocities and projected positions of tracer populations, such as dwarf galaxy satellites, 
globular clusters, distant halo stars and planetary nebulae or PNe (see e.g., Wilkinson et 
al. 2002 for a recent review). These do not follow the underlying mass distribution which 
is dominated by the dark matter. For nearby galaxies, the advent of wide-field imaging and 
efficient planetary nebulae spectrographs (e.g., Mendez et al. 2001; Halliday et al. 2002) 
means that the quality and quantity of such datasets are destined to improve dramatically 
in the next few years. 

The purpose of this paper is to present a new and simple mass estimator for galaxy 
haloes tailored for tracer populations and to apply them to new datasets for the Andromeda 
(M31) galaxy. It is the generalisation of the projected mass estimator to the commonest 
case of all - namely that when the number density of the tracer population is different from 
the overall mass density. 



On dimensional grounds, any mass estimator derived from projected data on a set of N 
tracer objects involves weighted means of viosf-Ri/G, where viosi is the line of sight velocity 
(relative to the mean or systemic velocity) and Ri the projected position of the ith object 
relative to the center of the galaxy. The projected mass estimator is 



This was introduced by Bahcall & Tremaine (1981) for test particles orbiting point masses 
(such as stars around a black hole or distant companion galaxies orbiting a central massive 
galaxy). The constant Cproj = IG/tt for test particles with an isotropic velocity distribution 
orbiting a point mass, whereas Cproj = 32 /tt for test particles moving on radial orbits. 
Subsequently, Heisler et al. (1985) extended the projected mass estimator to the case in 
which the galaxy members are believed to track the total mass, such as galaxy groups. 
The constant Cproj = 32/7r for particles with an isotropic velocity distribution, whereas 
Cproj = 64/7r for particles moving on radial orbits. 

We wish to generalise the mass estimator (1) to the case when the objects are drawn 
from a tracer population (such as globular clusters or PNe) between radii Tin and rout- Let us 
assume that the tracer population is spherically symmetric and has a number density which 
falls off like a power-law 
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within the region of interest. The underlying gravity field is assumed to be scale-free 



^(r) = { (3) 
—Vq log r if q; = 0. 

When a — 1, this corresponds to test particles orbiting a point-mass; when a — 0, the 
satellites are moving in a large-scale mass distribution with a fiat rotation curve; when 
a = 7 — 2, the satellites follow the total self-gravitating mass (the case studied by Hcislcr ct 
al. 1985). Properties of such tracer populations in scale-free spherical potentials have been 
deduced by Evans, Hafner & de Zeeuw (1997). 

Suppose our tracers are distributed over a range of projected radii i?in and i?out (cor- 
responding to the three-dimensional radii r^^ and rout)- What mass distribution do they 
probe? By Newton's theorem, there is no knowledge whatsoever on the mass distribution 
outside Tout- Accordingly, the best we can do is to estimate the total mass within the radius 
of the most distant satellite, namely 

M = '-^rlX, (4) 

from the dataset of positions and velocities. Let (r, 9, 0) be standard spherical polars with 
the z axis defining the line of sight, so that the projected position R is rsin^. Denoting the 
phase space distribution function of the tracer population by /, then the average value of 
vios^R in the survey is 

(vios'^R) -Jljj Mos'^Rdvdr = — y drj^ dO pa^r^ sin^ 9. (5) 

Here, we have computed the average over the ensemble between the radii Tin and rout probed 
by our data. In this formula, o"ios is the line of sight velocity dispersion, while Mx is the mass 
in the tracer population, which is given by 

Mt = <( (6) 
47raVo log(rout/nn) if 7 = 3. 



2.1. Isotropic Populations 



For the moment, let us assume that the velocity distribution of the tracers is isotropic 
so that the radial velocity dispersion cr^ is the same as the tangential Ut- In this case, Evans 
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et al. (1997) solve the Jeans equations and derive the result 



For example, if the potential is isothermal (a = 0) with a fiat rotation curve of amplitude Vo 
and the density of the population falls off like r~'^, then the velocity dispersion is just Vq/ ^ 
- which is a well-known result from the last century (e.g., Smart 1938). Now, substituting 
(6) and (7) into (5), we obtain: 

, , _ ^Vo^g ,3 V""" - ™^ 3-7 l-(rin/ro.t)^— ^ 

- Mt(«+7) X, \r) "4(«+7)4-«-7 l-(nn/rout)^-^ ■ ^ ^ 

In other words, suppose we gather positions and velocities for a population of globular 
clusters or PNe with a three-dimensional number density falling like between an inner 
radius and an outer radius rout- The mass enclosed within the outermost datapoint is 
given by the simple formula 

i 

with 

^_ 4(a+7)4-a-7 l-(/-iu/^ouJ'~^ 

TT 3-7 l-(rin/rout)^' ^ ' 

We refer to this as the tracer mass estimator. At first sight, it seems that the constant C 
in eq (10) is singular when the tracer population density falls off like (that is, 7 = 3) 
or when a + 7 = 4 (for example, if the tracer population density falls off like in an 
isothermal (a = 0) potential). However, a careful application of L'Hopital's rule shows that 
the constants are well-defined. They are listed in Table 1 for convenience. 

How do we set the parameters occurring in eq (10)? It is always reasonable to assume 
that Till ^ where Ri^ is the projected radius of the innermost datapoint. If the dataset 
is derived from a wide-angle survey, in which the population is traced out to large radii, 
then it is also reasonable to assume that rout ~ -Rout, the projected radius of the outermost 
datapoint. But, this assumption is not appropriate for a dataset restricted to the inner parts 
only, because this will usually contain objects at larger three-dimensional radii projected 
into the sample. The parameter 7 can be calculated from the surface density of the tracer 
population between Rin and it!out- However, a (and also possibly rout) need to be set using 
astrophysical considerations. We are envisaging applications to the outer parts of galaxies 
probed by globular clusters and PNe, and so it is reasonable to set a ~ as the potential 
is probably close to isothermal. For example, Wilkinson & Evans (1999) studied the Milky 
Way galaxy and found that a near- isothermal potential is valid out to ~ 170 kpc. 
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Table 1: The constant C in the tracer mass estimator eq (9) for the two special cases (7 = 3 
and a + 7 = 4). The number densities of tracer populations belonging to the spheroid or 
stellar halo are often oc or oc in the outer reaches, so the special cases often occur in 
practice. 



Case Constant 



^, _o ^ _ 4(q;+3)(1-q;) log(rout/r„.^ 
7 — o o — ^ 1 /'^ ^ i-a 



l-(rin/rout) 

3-7 



7r(3-7) log(rout/nn) 



Comparison of the result (9-10) with previous work is instructive. For example, when 
q; = 1, the potential is Keplerian and we find that 

i 

When q; = 7 — 2, the tracer population follows the underlying dark matter density and we 
find that 

These are not quite the same as the results given by Bahcall & Tremaine (1981) and Heisler 
et al. (1985). These investigators assumed that the sample was gathered from the center 
of the galaxy or cluster to infinity. In a number of integration by parts, boundary terms 
could therefore be legitimately dropped under the assumption pal — as r — > and as 
r — > oo. In our calculation, we have always performed the averaging over a finite range of 
radii and the contribution from the boundaries does not generally vanish. This is the origin 
of the difference. Even when we formally take the limit Tin and rout oo, the boundary 
terms still do not vanish as the density distribution (2) of the tracer population is cusped. 
Of course, such a limit is purely formal, as our estimator is derived under the assumption 
that the power-law density distribution holds over a certain regime only. 



2.2. Anisotropic Populations 



Now let us extend the calculation to deal with populations with anisotropic velocity 
distributions, although still retaining the assumption of spherical symmetry for the density 
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distribution. Binney's (1981) anisotropy parameter p — 1 — af/al is often used to mea- 
sure the relative importance of the radial and tangential ctj velocity dispersions. When 
/3 — > — cxD, this is the circular orbit model and when /3 = 1, this is the radial orbit model. 
For models in which the anisotropy (3 does not vary with radius, the line of sight velocity 
dispersion a\os is related to the radial velocity dispersion by 



aL = a,^(l-/3sin2^). (13) 



Substituting into (5), we find 



^"^-'^^ = ^^-lMp£^^ (14) 
The radial velocity dispersion cr^ of the tracer population is (Evans et al. 1997) 

Prom this, we deduce 

M = 

GN 

with 



^_ 16(a+7-2/3)4-a-7 l-{r-Jr,^,f-'y 

7r(4-3/?) 3-7 l-(rin/rout)^-"-^ ^ ^ 

This is the extension of the result to populations with constant anisotropy. Unless there 
are compelling reasons to the contrary (e.g., a highly flattened system supported by an 
anisotropic velocity dispersion), we advocate assuming isotropy. The value of eq (16) is 
that it enables us to calculate the error in making such an approximation. For example, 
suppose we consider tracer populations with a number density falling like in an isothermal 
potential. The anisotropy of stellar populations is not usually more extreme than 2:1. 
Therefore, cq (16) tells us that mistakenly assuming isotropy leads to mass underestimates 
(overestimates) in the case of radial (tangential) anisotropy of ~ 30%. 



2.3. Monte Ccirlo Simulations 

We test the performance of the mass estimator with Monte Carlo realisations of exact 
solutions of the collisionless Boltzmann equation. The phase space distribution functions 
corresponding to power-law density profiles in power-law potentials are given in Evans et al. 



- 7- 



Tablc 2: The percentage of mass estimates in error by a large factor (too big or too small 
by 50 per cent). The numbers are derived from 10000 Monte Carlo simulations of isotropic 
tracer populations with a number density falling like in an isothermal potential. The 
sample size N and extent rout/^n are chosen so as to be typical of datasets of satellite galaxies 
(upper panel) and globular clusters (lower panel). 



Estimator 


Number 


^out/'^in 


Too Low 


Too High 


Tracer 


10 


100 


20% 


37% 


Projected 


10 


100 


53% 


11% 


Virial 


10 


100 


87% 


0% 


Tracer 


100 


10 


2% 


2% 


Projected 


100 


10 


72% 


0% 


Virial 


100 


10 


100% 


0% 



(1997). To generate Monte Carlo realisations, the speed is picked from the distributions 

r ^2-2, |^(,) _ l^2|P7-3a-2,(2-a)]/(2a) ^ ^ 
f{v) OC I ^ (18) 

I t;2-2/3exp(-^) ifa^O. 

For a > 0, the maximum velocity at any position is •\/2'0(r); for a < 0, the velocities can 
become arbitrarily large. Following Binney & Tremaine (1987), let us introduce spherical 
polar coordinates in velocity space (t>,^,?7), so that the velocities resolved in spherical polar 
coordinates with respect to the center are then 

Vr = vcos'rj, Vd = V sin rj cos ^, ~ vsmrjsinC,. (19) 

If ^ is generated uniformly in [0, 27r], while rj is picked in [0, tt] from the distribution 

F(r/) OC |sin?7|i+2^ (20) 

then the velocities are generated with the correct anisotropy. Finally, a random projection 
direction is chosen and the line of sight velocity vios and projected position R calculated. 
This gives us synthetic datasets of anisotropic spherical tracer populations in spherical haloes 
with which to test the performance of our mass estimators. 



Figure 1 shows the medians and quartilcs for the mass estimates of 10000 samples 
constructed from Monte Carlo realisations. We test the tracer mass estimator for both 
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isotropic (eq. 9) and anisotropic (eq. 16) populations, as well as the projected and virial 
mass estimators in the forms 

For each of the mass estimators, four models are shown in each panel. In the upper two 
panels, the tracer population has a number density falling like r^^ and is generated between 
i?in and i?out = 10i?in, while the rotation curve is fiat (a = 0), The models vary in the 
anisotropy of the velocity distribution (either ctj : cr^ = 1 : 1 or (Jt : cr^ = 2 : 1 or cr^ : (7^ = 1 : 
2). In upper panel, the number of objects is N — 10, which is typical for a sample of sateUite 
galaxies. In the middle panel, the number of objects is = 100, which is reahstic for a 
sample of globular clusters with radial velocities. Finally, in the lower panel, the rotation 
curve is falling in a "half-way Keplerian" manner (a = 1/2) and the tracer population is 
falling like r^^. The sample size is again N = 10. The tracer mass estimator outperforms 
the projected and virial mass estimators in every case. This is as it should be, since the 
tracer mass estimator has been devised for the precise purpose of mass estimation from 
tracer populations, while the other two estimators have been devised with other applications 
in mind. It is interesting to ask how often the estimates are too low or too high by 50%. 
This information is recorded in Table 2 for some of the simulations. The projected and virial 
mass estimators yield systematic underestimates if they are mistakenly applied to tracer 
populations. Further, the situation is not improved by increasing the size of the tracer 
population - the estimators merely converge to an underestimated mass. 

Finally, we test the tracer mass estimator against synthetic datasets drawn from a self- 
consistent Plummer model with a core radius b. This is a important thing to do as the tracer 
estimator has been derived for, and tested against, the scale-free case. The Plummer model 
is only approximately scale-free well outside the core region (r ^ 6). The potential, density 
and distribution of speeds are 

P"^ + (^2+^^2)1/2 ' f{v)(xv\ij{r)-^vy/^ (22) 

Figure 2 shows how the estimators fare for sample sizes N — 10 and N — 100 using datasets 
gathered in the outer parts (it!in = 106, it!out = 1006) and including some core contamination 
{Rin — b,Rout = 1006). For a self-consistent population, the projected mass estimator 
performs well, as it should do. Nonetheless, the tracer mass estimator (with 7 = 5, « =1 
and j3 = 0) gives slightly better results. The reason why the tracer mass estimator still 
out-performs the projected mass estimator is that the former includes a correction for the 
limited radial coverage of the sample, whereas the latter does not. 
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3. Application: the Andromeda Galaxy 



It is important to know the mass of the Andromeda Galaxy (M31) for a number of 
reasons. First, the density profiles of dark matter haloes are an important constraint on 
cosmological theories of galaxy formation (e.g., Sellwood 2001). The haloes of nearby galaxies 
are particularly important as they can be studied in much greater detail. Second, the current 
generation of pixel lensing experiments towards M31 (e.g., Auriere et al. 2001, Pauhn- 
Henriksson et al. 2002) requires an accurate estimate of the total mass so as to constrain the 
fraction in compact objects capable of producing microlensing events. We shall consider two 
datasets which probe Andromeda's dark halo, namely the globular clusters and the dwarf 
spheroidal satellite galaxies. 



Perrett et al. (2002) give the positions, velocities and metaUicities for over 200 globular 
clusters in M31 out to a projected distance ~ 30 kpc. The data were obtained using fiber 
optic spectroscopy on the William Herschel Telescope and the typical errors on the radial 
velocities are ±12 kms~^. When combined with earlier data from Huchra, Brodie & Kent 
(1991), this gives a grand total of ~ 300 globular clusters. From the bimodality of the metal- 
licity distribution, Perrett et al. show that the M31 globular clusters fall into metal-poor 
and metal-rich categories with different kinematics. As judged from the double Gaussian 
fit to the metaUicity distribution, there are 177 globular clusters with a probabihty of be- 
longing to the halo of greater than 90%. The cumulative number distribution of the halo 
globular clusters is shown in Figure 3. Only beyond 30 arcmin does the projected number 
density of the halo globulars fall off like a power-law, namely R^^. Accordingly, we work 
with the 89 halo globular clusters with projected radii greater than 30 arcmin. The mean 
rotation amplitude of this sub-sample is (v^) ^110 kms^^. Adopting a distance of 770 kpc, 
the globular clusters lie in projection between 6.8 kpc and 33 kpc. However, the globular 
cluster population in M31 certainly extends out to ~ 100 kpc (e.g., Hodge 1992). Using 
Monte Garlo simulations to build samples of 88 globular clusters with 6.8 kpc < i? < 33 kpc, 
we find that typically 20% of the sample have r > 33 kpc with the outermost datapoints 
typically at r ~ 100 kpc. Accordingly, in the formula for the tracer mass estimator, we set 
a — (isothermal-hke galaxy), 7 ~ 4, rin = 6.8 kpc and rout ~ 100 kpc. 

The Jeans equation for a population with constant anisotropy (3 about a mean velocity 
(v^) is (e.g., Binney & Tremaine 1987, section 4.2) 



3.1. 



Globular Clusters 




(23) 
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The first term on the right-hand side describes the contribution of rotation to M(r), the 
second term the contribution of pressure. With obvious notation, we write 



rot / \ 2 ^-'"press 2 



5 log pa; 
dlogr 



+ 2/3 



(24) 



r r 

We compute Mrot ~ 3 x IO^-'^Mq directly from the mean rotation amplitude of the sample. 
We use the tracer mass estimator, applied to the observational data with the mean rotation 
velocity subtracted from the line of sight velocity, to compute Mprcss ~ 9 x IO^^Mq assuming 
isotropy. Hence, the total mass within ~ 100 kpc of the center of M31 is ~ 1.2x 1O^^M0. Note 
the tracer mass estimator applies to a pressure-supported tracer population and therefore 
any net rotation of the system must be subtracted before the velocities are used to estimate 

-^press • 

How docs the result from the tracer mass estimator compare with earlier work? Evans 
& Wilkinson (2000) give the results of more sophisticated modelling. Here, the phase space 
distribution function of the globular cluster and the satellite galaxy population is built 
and then convolved with the errors to find the probability of the dataset given the model 
parameters. This is inverted using Bayesian likelihood techniques to give the probability of 
the enclosed mass given the positions and velocities of the halo globular clusters. Using the 
combined globular cluster and satellite galaxy sample, they reckoned that the most likely 
total mass of M31 is ~ 12^g^ x IO^^Mq. An independent constraint on the mass of M31 
within 30 kpc is given by the HI rotation curve. This was measured to distances of ~ 30 kpc 
along the major axis by Newton & Emerson (1977). Their estimate of the asymptotic circular 
speed is 230 kms~^. Assuming the circular velocity curve maintains this amplitude out to 
100 kpc, then the mass within 100 kpc is 1.2 x lO^^M©. We conclude that the results are all 
consistent and support a picture in which M31's dark halo extends out in an isothermal-like 
manner to at least 100 kpc. 



3.2. Satellite Galeixies 

There are 15 probable companion galaxies to M31. Their galactic coordinates, projected 
position from the center of M31 and line of velocity are listed in Table 3. The observed radial 
velocities contain contributions from the Galactic rotation and from the relative velocity 
between the Galaxy and M31. Assuming that the transverse velocity of M31 with respect 
to the Galaxy is zero, we remove these contributions and list the corrected velocities in the 
frame of M31 in Table 3. 

The number density of M31's sateUite galaxies falls hke r~^-^ (see e.g., Evans & Wilkinson 
2000). Notice that the projected positions of many of the objects lies well beyond the regime 
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in which the M31 halo is expected to be isothermal. Rather, the underlying gravitational 
potential is probably dominated by the monopole component at such huge distances. This 
suggests that it is appropriate to apply the tracer mass estimate with a ^ 1 and 7 ^ 3.5 
to the satellite galaxy datasct. This gives the mass of M31 as ~ 1.1 x IO^^Mq. This is in 
good agreement with the more sophisticated modelling in Evans et al. (2000), who claimed 
that the total mass is 7 — 10 x lO^^M© from the same dataset. In fact, it is in excellent 
accord with still earlier work by Bahcall & Tremaine (1981), who also found a mass of 
~ 1.0 X 10^2^0 from just five companions (M32, M33, NGC 147, NGC 185 and NGC 205). 
Despite the seeming agreement, the uncertainty in the mass is at least a factor of 3. For 
example, if instead we choose a ~ 0.5 (suitable for a halo with a rotation curve slowly falling 
hke oc r~°^^), then the mass of M31 as deduced from the tracer estimator more than triples 
to ~ 3.8 X IO^^Mq. 

An assumption underlying the tracer mass estimator (as well as the projected and virial 
estimators) is that of a steady state equilibrium. This is reasonable enough for globular clus- 
ter and PNe datasets, but probably not for satellite galaxies. The dSph and dirr companions 
of M31 may well be falling in for the first time, in which case the assumption of a steady 
state is probably not a valid description of the dynamics. It would be valuable to calibrate 
the importance of the effects of disequilibrium against high resolution N-body simulations 
of the Local Group. 

4. Conclusions 

The main accomplishment of this paper is the introduction of the tracer mass estimator. 
This is a new and simple way of estimating the enclosed mass from the radial velocities and 
projected positions of a set of objects. It is tailored for the case of tracer populations such as 
globular clusters and halo stars whose number density falls off hke a power-law with distance 
(at least within the range covered by the survey). The estimator works by averaging the 
quantity (constant) x (projected distance) x (radial velocity)^/G over the sample. The value 
of the constant is straightforward to calculate. It depends on the radial fall-off of the tracers, 
the inner and outer radii marking the range of the sample, the velocity anisotropy of the 
population and the underlying gravitational potential. In the absence of further information, 
we recommend assuming isotropy and isothcrmality for estimating galaxy masses, at least 
in the case of globular cluster and planetary nebulae (PNe) datasets. The inner and outer 
radii can usually be inferred from the sample itself, sometimes using additional astrophysical 
evidence. 

We have applied the tracer mass estimator to the globular cluster and satellite galaxy 
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Table 3: Data on the companion galaxies of M31 taken from Grebel (2000) and Evans et 
al. (2000). Listed are Galactic coordinates (£,6), the projected distance from the center of 
M31 R in kpc, corrected line of sight velocities vios (adjusted for the solar motion within the 
Galaxy and the radial motion towards M31) and object type. (Note: the coordinates for the 
CasdSph given in Table 1 of Evans et al. (2000) are incorrect). 



TVT 

Name 









R 


Vlos 


Type 


M31 


121.2 


-21.6 


— 


— 


Sbl-II 


M32 


121. 


,1 


-22.0 


5 


+95 


E2 


NGC 205 


120. 


,7 


-21.1 


8 


+58 


dSph 


NGC 147 


119. 


,8 


-14.3 


100 


+ 118 


dSDh/dE5 


NGC 185 


120. 


,8 


-14.5 


95 


+107 


dSph/dE3 


M33 


133. 


.6 


-31.5 


198 


+72 


ScII-III 


IC 10 


119. 


.0 


-3.3 


243 


-29 


dirr 


LGS3 


126. 


.8 


-40.9 


262 


-38 


dIrr/dSph 


Pegasus 


94. 


,8 


-43.5 


397 


+86 


dIrr/dSph 


IC 1613 


129. 


,7 


-60.6 


489 


-58 


IrrV 


Audi 


121. 


,7 


-24.9 


44 


-85 


dSph 


And II 


128. 


,9 


-29.2 


138 


+82 


dSph 


And III 


119. 


,3 


-26.2 


67 


-58 


dSph 


AndV 


126. 


.2 


-15.1 


07 


-107 


dSph 


And VI 


106. 


.0 


-36.3 


260 


-65 


dSph 


Gas dSph 


109. 


,5 


-9.9 


215 


+23 


dSph 



datasets of the Andromeda galaxy (M31). The estimates suggest a picture in which MSl's 
dark halo is isothermal out to at least ~ 100 kpc and has a mass of ~ 1.2 x IO^^Mq. 
These estimates are in good agreement with the results of more sophisticated modelling. 
The advantage of the tracer mass estimator is that it is so much quicker to apply and 
requires much less effort. It should be used in preference to the projected mass estimator 
for populations in the outer parts of galaxies. One of the fundamental assumptions of the 
projected mass estimator (at least in the form which it is usually encountered) is that the 
population shadows the overall mass density. This assumption breaks down for populations 
in the outer parts of galaxies. We envisage the tracer mass estimator as being particularly 
useful for early-type galaxies for which there is no gas rotation curve and so mass estimates 
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almost always derive from dynamical modelling of tracer populations. 
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Fig. 1. — Medians and upper and lower quartiles for the mass estimates of 10000 tracer 
samples constructed from realisations of exact scale-free models. Upper Panel: The tracer 
population has a number density falling like in an isothermal potential. The data is 
gathered between Rout — 10i?in. Prom left to right, the models are: (i) 10 objects drawn 
from an isotropic velocity distribution, (ii) 10 objects drawn with at : ar — 2 : 1, and (iii) 10 
objects with at : — 1 2. For anisotropic populations, knowledge of the anisotropy P is 
assumed and the tracer mass estimator in the general form (16) is used. Middle Panel: As 
the Upper Panel, but the number of objects is now 100. Lower Panel: As the Upper Panel, 
but the tracer population has a number density falling like in a potential with a falling 
rotation curve ip oc r~^/^. 
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Fig. 2. — Medians and upper and lower quartiles for the mass estimates of 10000 tracer 
samples constructed from realisations of exact Plummer models with core radius b. The 
self-consistent population has a number density falling like in a potential falling like r"^ 
at large radii. Again, the performance of the tracer, projected and virial mass estimators is 
compared. From left to right, the models arc: (i) 10 objects gathered between i?in = 106 
and _Rout = 1006, (ii) 10 objects between = b and -Rout = 1006, (iii) the same as (i) but 
with = 100, and (iv) the same as (ii) but with = 100. 
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Fig. 3. — Cumulative number of M31 halo globular clusters versus projected radius R (in 
arcmin). The dotted line shows the profile expected for a population with three dimensional 
density distribution p ~ r~** between 30' and 100'. 
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